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ABSTRACT 


A perturbation guidance scheme is developed to keep a 
tilt-rotor VTOL aircraft close to a predetermined nominal 
flight path during take-off and landing. 

A simulation of the guidance scheme applied to the 
Bell Model 266 tilt-rotor VTOL gave satisfactory behavior 
in the presence of initial errors and wind disturbances. 
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CHAPTER I 
Introduction 

Vertical take-off and landing (VTOL) aircraft offer a 
potentially useful means of short-haul transportation. 

Since they do not require long runways for take-off and 
landing, it is possible to provide convenient air service 
using only small land areas for terminals. 

Several different kinds of VTOL aircraft have been built 
and flown in the past decade. The tilting proprotor VTOL 
aircraft promises to be among the most efficient and least 
noisy during transition and hover. However, >t has only 
a moderate cruise speed. 

For economic viability, VTOL aircraft must be capable 
of operating (a) at night, (b) in fog, haze and cloudy 
weather, (c) near tall buildings, and (d) in a crowded air 
traffic environment. Because of the low departure and ap- 
proach speeds of these aircraft there may be many aircraft 
in the vicinity of the airport during peak travel periods. 
Therefore, it is important to devise a good guidance scheme 
for these aircraft in the terminal area. 

The present report describes a guidance scheme for a 
tilt-rotor aircraft during take-off and landing. The scheme 
has been simplified as much as possible to permit its 
implementation with an onboard computer. The Bell Model 266 
is taken as the example aircraft. 

Chapter II develops a mathematical model for the 
longitudinal motion of a tilting proprotor VTOL aircraft. 

Chapter III describes the technique used for finding 
guidance laws during take-off and landing. 
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B. CRUISE MODE 

Figure 1.1 The Model 266 Tilt-Rotor VTQL Transport 
Designed by Bell Aircraft 
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Guidance laws are given for take-off in Chapter IV and 
for landing in Chapter V. The guided aircraft behavior in 
the presence of initial and wind disturbances is studied in 
both chapters. 

Chapter VI develops a mathematical model for the air- 
craft lateral motions. A guidance law is developed and 
guided aircraft behavior is studied in the presence of 
lateral winds. 

Appendix A gives relevant data for the Bell Model 266. 

Appendix B summarizes the longitudinal equations of 
motion and the weighting matrices used for computing guidance 
laws in Chapters IV and V. 
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CHAPTER II 


Longitudinal Equations of Motion for 
a Tilt-Rotor VTOL Aircraft 


2,1 Introduction 

In order to develop guidance laws for take-off, landing 
and other flight conditions, it is necessary to develop a 
mathematical model, which is accurate enough for this pur- 
pose. Modeling aerodynamic forces is complicated by (a) the 
interaction between the rotor downwash and the wing, and 
(b) the change in configuration with rotor tilt. 

The model developed here is fairly general and can be 
used for most tilt-rotor VTOL aircraft. The Bell Model 266 
is taken as an example in our discussion. The Model 266 is 

a high-wing, twin engine, tilting proprotor aircraft. It 

_ 5 

has a disc loading of 58.6 Kg m and a wing loading of 358 
_ 2 

Kg m at its design gross weight of 12,700 kilograms. 

The two nacelles are located at the wing tips and can 
be rotated through slightly more than 90 degrees. Each 
nacelle accomodates an engine, transmission and other ac- 
cessories necessary to drive a proprotor. The proprotors are 
similar to the rotors in a tandem rotor helicopter. Each 
rotor has collective pitch control and cyclic pitch control 
in the longitudinal and the lateral directions. The two 
engines are connected through an intermediate transmission 
so that in case of one engine failure, the two rotors can 
be driven at a reduced power. 

A compromise is made between Tow disc loading for ef- 
ficient hover and a high disc loading (in the propeller 
mode) for efficient cruise. A high wing loading is chosen 
to obtain good performance at cruise speed. 
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Some parameters of the Bell Model 266 are given in 
Appendix A. 

2.2 Aerodynamics of the Rotor 

We first develop expressions for the thrust and in-plane 
forces acting on the rotor and the rotor downwash under 
different flight conditions. 

2.2.1 Thrust 

The shaft horsepower required to generate a certain 
thrust under given flight conditions depends on the blade 
design. Different blade twists are required for efficient 
hover performance than for efficient cruise performance. 

In our analysis, we use a blade twist which is efficient in 
cruise for two reasons: (a) a twist distribution which is 
efficient for cruise suffers a smaller penalty in hover than 
vice-versa, and (b) the aircraft is in the cruise mode much 
longer than it is in the hover or transition mode. 

The thrust available and power required at various col- 
lective pitch settings and flight conditions are obtained by 
using a mixture of Blade Element Theory and Momentum Theory 
As a first approximation, the angle between the free stream 
velocity and the normal to the tip path plane is assumed 
zero. Later analysis shows that when the free stream 
velocity is substantial this inflow angle is quite small. 

A correction can be applied for finite values of this angle. 
Fig. 2.1 shows a velocity diagram for a section of the blade 
at radius r. If R is the rotor radius and Vy the tip speed 
of the blades, the speed of this element of the blade is 
jj- Vy. Let V be the free stream velocity and v x the induced 
velocity at radius r. v is parallel to the lift on this 

A 

element. The rotor downwash is supposed to be free of ro- 
tation. Let y be the twist at radius r and $ the col- 

X V 
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Figure 2.1 Velocity Diagram for Rotor Blade 
Element at Radius r 

lective pitch. Then the angle of attack, a x , on this blade 
element is given by 

“x = *c +Y x" Y i" Y » (2.2.1) 

where y is the inflow angle and y^ is the induced inflow 
angle. From Figure 2.1 it is clear that 
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Y.j = arcsm 


((fV T )SV^) 


The resultant wind speed V R , at this blade element is 
given by 


Vr * 


( 2.2.4) 


If c is the blade chord (assumed constant), the lift, 
dL on this element of the blade can be determined 


dL = P C L (« X )V R cdr 


7 P ^L {a x ] ^ + V^ -v x j cclr 


(2.2.5) 


Similarly the drag, dD on this blade element is 


dD * 7 0 C D x '«x) V R cdr 


7 p (« x ) ^ j) J cdr (2.2.6) 


If n is the number of blades, the total lift and drag 
produced by the section of rotor disc between radius r and 
r+dr can be expressed as 


\ p ncR C L (« x ) ^(rVj) 2 +V 2 -v x j d(r) 


(2.2.7) 
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dD d = \ P ncR Cp (a x ) 
X 


(rV T ) 2 *V 2 -v2] 


d(r) 


( 2 . 2 . 8 ) 


Resolving these forces perpendicular and parallel to 
the rotor disc gives the thrust component, dT and the torque 
component, dD T , respectively: 


dT = dL D 

COS 

(y+y^) - dDp 

sin 

(y+y-j ) 

(2.2.9) 

dD T = dLp 

si n 

(y+y^ ) + dDp 

cos 

(Y + Y-j ) 

(2.2.10) 

be the 

tip 

loss factor. 

Then 

the total 

thrust 


on N rotors is given by 


r BR 

T = N / dT 

J A 


J r \ p ncR N j(rVy) 2 +V 2 -vj|j C L (a x ) cos ( y+y ^ ) 

■'] 


-C n (a ) sin ( y+y -5 ) I d ( r ) 
u x x 


Ip ncR V* N /V + 4 - %) 

c 1 J 0 V Vy Vy ’ 


C (a x ) COS (y+Y^) “ Cp (o x )sin(Y+Y 1 * ) 


dr 


(2.1.11) 


Gessow and Myers (GE-1) give a convenient expression for 
the tip loss factor in terms of thrust coefficient, Cy 
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and number of blades n 


B = 1 - 




n 


( 2 . 2 . 12 ) 


The torque produced by this section of the disc can be 
computed 

dx = rdD 

t 

The total torque on one rotor is 


T 



( 2 . 2 . 13 ) 


Hence the power, Q, required to drive the N rotors at 
this collective pitch setting becomes 


Q = xfi N 




rf! y p N n c R 


; 2 v 2 +v 2 2 

r V T +V -v x 


Cl (a x ) sin (y +Y 1 > ) 


+ C D (a x ) cos (y+y . ) 


dr , 


( 2 . 2 . 14 ) 


where is the rotational speed. Also, 


rft = •£ Rft = rVj 
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« " /’ 7 PNncR V^V^-v 2 ] 


• |c^ (a x ) sin (Y + Y-j) + Cp (« x ) cos (v + Y^)Jdr 


= j pNncR V 




A ^ 

v ? v r 


|c L (c» x ) sin (y+y^)+C d (a x ) cos (Y + Y^dr 

v. X X 


(2.2.15) 


Let p be the average pressure difference across the 

A 

disc at radius r. Then, 


C. & ncR (r 2 v!+V 2 -v 2 )dr 

- dL L x _ 1 x 

x ‘ d( Area ) " 2irr dr 


= i P C. (r 2 V 2 +V 2 -v 2 ) (2.2.16) 

c L x 2 ttR> 1 x 


We can find v by expressing p in another way. By 

A A 

momentum theory 

p x = l (V r v2) * (2.2.17) 

where V f is the maximum airspeed in the rotor downwash. Since 
the maximum induced velocity is twice the induced velocity 
at the disc [GE-1], V f can be expressed as 
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Vf . V + 2* x 


or. 


V 2 * V 2 +4v 2 +4Vv x cos (y + Y-j) 


(2.2.18) 


Using equations (2.2.16) to (2.2.18), we have 


I 4v x +4Vv x cos (y+y^ 


1 

7 p 


ncR tlv 2^1/2 2\ 

I rr(rV T +V -v ) 

L x 2irR*r 1 x 


ncR 



a, solidity ratio 


(2.2.19) 


v 2 

X 


1 + 


C. CT 
x 

8r 


+ v x V cos (y+y^) 


[r 2 vf+V 2 ] = 0 
8r 

( 2 . 2 . 20 ) 


This equation can be solved for v x at each radius. 

The following expressions were used for the lift coeffi 
cient and the drag coefficient on any blade section. 


C 


6.0 a x 



1 

T 

( 2 . 2 . 21 ) 

x not allowed 
(stall condition) 


C D X = .008 


.008 C, + .020 C 2 
L x L x 


( 2 . 2 . 22 ) 


2.2,2 Computation of In-Plane Forces 

In almost all flight conditions considered here the 
component of in-plane force in the forward direction is very 
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small. This force has a significant effect only during high 
speed flight in the helicopter mode. Thus, for purposes of 
this analysis, it will be sufficient to determine this force 
approximately. It is computed by taking an average drag 
coefficient Cpg on the blade. The drag on an element of the 
blade between radius r and r+dr situated at azimuth angle \ i> 

from down wind position is (Figure (2.2)): 

dD ” 7 V T +V Sln * cos cdr c d (2.2.23) 

B 

Its component in the downwind direction is, 

dH = 7 p (fr V T +V sin * cos cdr c d sin * (2.2.24) 

B 

Average total drag on all blades can now be determined. 



Figure 2.2 In-Plane Forces on a Blade Element 
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i /* 2ir /*R i 

H ' ™o Jo C|) B 5ln * 


■(* h 


+ V sin ib cos 


>) 2 <4 


A 

/•l , / 2rV T V cos 4 \ a 

7 . 7 p ncR c D r \ 2 ) dr 


, V T V cos 4 > 

7 p ncR C D g 2 


Hence drag force on N rotors is, 

H = ■J- p o NitR 2 C d (V t V cos 

B 


(2.2.25) 


(2.2.26) 


This is an approximate formula and corresponds to the 
first term in equation (4.38) of Nikolsky ( N I - 1 ) 


2.2.3 Induced Velocity 

The computation of wing drag and lift requires know- 
ledge of rotor downwash. The downwash velocity depends on 
the radial distance from the axis of rotation and distance 
perpendicular to the tip path plane. To simplify calculation 
of lift and drag over the wing in the next section, the 
induced velocity over that part of the wing immersed in 
rotor downwash is assumed to be constant and equal to the 
maximum induced velocity. Using momentum theory this 
velocity is 2v where v is obtained by solving the equation 
(6E-1), 


v[V 2 +v 2 +2Vv cos ii] 1 ^ 2 = — 


r 


N 2nR p 


(2.2.27) 
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2.2.4 Thrust for the Example Aircraft 

A computer program was written to determine the maxi- 
mum available thrust at rated power. Figure 2.3 shows the 
variation of thrust with inflow velocity. For nonzero 


Rotor 

Thrust 

[Kilo-Newtons] 



0 40 80 120 160 

Free Stream Air Speed [m sec’^ ] 


Figure 2.3 Variation of Rotor Thrust with 
Free Stream Air Speed 
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values of Inflow angle, the component of velocity in the 
direction normal to the tip path plane is used to read thrust 
from Figure 2.3. This gives conservative results. 

At least one author (HA-1) has suggested using trailing 
edge flaps on the blades to obtain almost optimal twist 
distribution along the blade under all flight conditions. 

The maximum available thrust for a rotor with blades having 
trailing edge flaps is shown in Figure 2.3 (dotted line). 

The maximum increase is about 4.5% over a rotor with flap- 
less blades. Notice that both perform equally well at 
cruise. The rotor with blade flaps performs better during 
hover. Such a minor improvement in performance with flaps 
does not seem to warrant the additional mechanical complica- 
tion. 

2.3 Equations of Motion in Forward and Vertical Directions 

The forces acting on a tilt rotor VTOL aircraft can 
be divided into two classes: (a) gravitational and (b) 
aerodynamic. The aerodynamic forces can be subdivided 
into two classes: (i) forces with all control elements in 
equilibrium position and (ii) additional forces generated 
by deflecting control surfaces. Gravitational forces are 
quite easy to find. The aerodynamic forces are approxi- 
mated by sums of forces on Individual components of the 
ai rcraf t . 

A simplified force diagram is shown in Figure 2.4. Let 
D be the drag and L the lift. The thrust and the in-plane 
force on the rotor are denoted by T and H respectively. The 
equations of motion parallel and perpendicular to the 
velocity of the vehicle mass center can be written as 

m ^ = -D-mg sin y + T cos ( n + 0 - y ) - H sin (n+e-y) 

(2.3.1) 
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mV = +L-mg cos y +T sin (n+0-y)+H cos (n+0-y) 

(2.3.2) 



Figure 2.4 Aircraft Reference Angles and Forces 

The drag and lift are expressed as sums of forces on 
individual components. 

D * + Dj + Dp + D|^q (2.3.3) 

L ■ L « + l t + l f + l nac (2 - 3 - 4) 

We consider these aircraft components one by one. 
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2.3.1 Wing 

The wing of a tilt-rotor VTOL airplane will usually 
have flaps and aileron-flaps. These flaps can be deflected 
through 60° or more during hover to reduce wing surface in 
the rotor downwash. The ailerons can also be used as flaps 
during hover and part of transition. 

A portion of the wing is in the rotor downwash. This 
has a significant effect during hover and the early part of 
transition where the induced velocity is high. It is assumed 
that the portion of the wing under the rotor disc is in the 
rotor downwash and faces a constant induced velocity. This 
induced velocity over the immersed part of the wing is 
approximated by the fully developed downstream velocity in 
the slipstream. Let e be the fraction of the wing in rotor 
downwash. This portion is denoted by subscript 1 and the 
other part by subscript 2. 

The angle of attack a over the portion of the wing 
in the free stream is given as 

o * 0 + a ^ - y (2.3.5) 

The flaps change the effective angle of attack. Let 
the effective angle of attack be o'. 

o' = a+ Aa^ ( 6f ) 

-a+g^6f (2.3.6) 

where <5f is the flap deflection. The lift coefficient 
depends on o' and the drag coefficient on a and 6f . 

°W 2 = 7 p C D w (a ’ 6f) v2 S W (l_e) (2.3.7) 

L w ■ \ P C L (o') V 2 S w (l- e ) (2.3.8) 
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Figure 2.5 Velocity Diagram over Immersed 
Part of Wing 

The downwash changes the magnitude and the direction of 
velocity. A simplified velocity diagram over the immersed 
portion of the wing is shown in Figure 2.5. The resultant 
velocity, V R , can be computed. 

V 2 = V 2 + 4v 2 + 4Vv cos 4 (2.3.9) 

The change in angle of attack, Aa, is obtained using 
the Sine Law of triangles. 



where , 


(2.3.10) 


hence, , ^ 

Aa = arcsin sin (<t>)J (2.3.11) 

The resultant angle of attack a R is written as, 


a R = a- Aa (2.3.12) 

Now the forces on the wing along the direction of V R , 
D^i , and perpendicular to the direction of V R , L^ R , are 
given by 1 

D w s j P C D ( a R , 6 f ) V R S w e (2.3.13) 

1 w 

L W B 7 P C L («'-Aa) vjj S w e (2.3.14) 

1 W 


These forces must be resolved along the direction of 

the free stream velocity and perpendicular to this direction 

to give D u and L u . 

W 1 W 1 

i i 

0,. = D,. cos (Act) + L u sin (Aa) 

W 1 W 1 W 1 

~ P V R e CCp ( o r » 4f) COS (Aa) 

U W 

+C, (a'-Ao) Sin Aa] (2.3.15) 

L W 
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L w ^ = -D^ sin (Ao) + cos (ia) 


“ "jr P Vp Sy e [*Cp (ct^tfi'f) si n( Aa) 

w 


+C, (a'-Aa) COS (Aa)] 

L W 


(2.3.16) 


The total drag on the wing can be written as. 


D w = d w 1 +d w 2 


L w " l w 1 +l w 2 


(2.3.17) 


Determination of Drag Coefficient and Lift Coefficient 

The flight regime of a tilt-rotor VTOL aircraft en- 
compasses angles of attack varying from -90° to 90° or more. 
These high angles of attack occur at low speeds, during 
take-off and during landing. Most of the high speed flight 
takes place at low angles of attack with most or all of the 
weight carried by the wing. A large error in drag co- 
efficient and lift coefficient at these high angles of at- 
tack will produce only a small error in the forces acting 
on the airplane because of low velocity. Therefore it is 
not necessary to have very accurate relations for and 
c L w for high positive or negative angles of attack. 

It is essential to find the positive and the negative 
angle of attack at which wing stall occurs. We suppose 
that with no flap deflection the stall occurs at angle of 
attack + a stall • The flap deflection increases the effec- 
tive angle of attack. It also increases the effective 
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angle of attack at which stall occurs. In terms of actual 
angle of attack the stall starts earlier because flap de- 
flection helps flow separation near the trailing edge. The 
actual positive and negative angles of attack at which stall 
occurs is approximated by using Figure 12 of reference 
(BE-1 ) 


+ 9f sf 

“s = “stall T 


(2.3.18) 


g f 6f 

°s = ’“stall' 


(2.3.19) 


The effective angle of attack a' is given by equation 
(2.3.12). In terms of this angle the stall condition can be 
written as 


. 9f Sf 

“s ' "stall + T~ 


(2.3.20) 


+ 9 f 6f 

“s “stall ~2 


(2.3.21 ) 


During the unstalled part of the flight a linear re- 
lation is assumed between the lift coefficient and the effec 
tive angle of attack, i.e., 

C l_(“ ' ) = a “' “s’ < a ‘ < “s + (2.3.22) 


As the angle of attack increase above a* the lift 
coefficient decreases. At a certain high angle of attack 
the lift coefficient becomes zero. This angle is approxi- 
mately 90° with no flap deflection but decreases with posi- 
tive flap deflection. It is reasonable to assume that the 
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lift coefficient is zero when the net projected area In 
the lateral direction is zero (Figure 2.6). 



Figure 2.6 Zero Lift Angle of Attack 


Thus , 


C L (a') = 0 


C - sin <$f 

tan { 7 ~ o) = c^C f -(l-cos 6 f ) 


and a similar expression for negative angles of attack 

C f sin 6 f 

C L U') = 0 tan(f -«) = - e ^r ^d r cos 6 f ) 
These equations are approximated for cr^ + 


C L ( o 1 ) = 0 


• -±f - f 


(2.3.23) 
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A linear relation is assumed between lift coefficient 
and angle of attack for a between stall angle and angle at 
which the lift coefficient is zero. Therefore, 


Cj^ ( o ' ) s aa ' 


i - i i + 

a s — a — a s 


aa 


,+ 


- £ 1 6 f + g f <5f -a' 


j ~ £—■ <5f + g^^f - a 


.+ 

s 


„ I % « 1 + 

a > a ^ 


aa ! 


C f 

^ + •£— Sf + a' g f 6f 


Tl , I 

7 * zr 


6f + a. 


W 




a 1 < a 1 

S 


(2.3.24) 


A plot of lift coefficient with angle of attack for 
no flap deflection and for one radian flap deflection is 
shown in Figure 2.7. 

In the unstalled regime the drag coefficient is given 
by 


C (a* ,6f) = C D («f) + C D Cf 
U W U W 0 U w 1 L 


I - I I + 

a < a < a 

S - _ s 


i 2.2 


* C n (6 f) + C n a'‘a 
\ W 1 


(2.3.25) 


The first term is the parasite drag which depends upon 
the flap deflection and the second term is the induced drag. 
A good approximation to the first term is (PE-1) 
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Figure 2.7 Lift Coefficient vs. Angle of 
Attack on Wing 
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D (6f) - 9 1 + 6 2 («f) 
W 


In the stalled regime the lift reduces; however the 
drag keeps increasing. The following expressions are used 
for drag coefficients for the stalled portion of the flight. 


C D (<x' .if) 


C D U (sf)+C D u C L + ®3< a ' -a s + ) 
W q W ^ max 


c n (6f)+c n C 2 +6,(<'-a') 
u w 0 U W 1 L max J s 


a 


I 


< a 


I ~ 

s 


The expression for Cp(a' ,6f) can be summed up as 



,6f) = B 1 +s 2 6f 2 +C D a 2 c^ 

W 1 


2 

+6 3 {o s’" a ' } 


a 


I 


< a 


i ” 

S 


6 1 + 6 2 6f 2 + C D a 2 ot ' 2 
W 1 


I “ I 

°s l a 


< a 1 + 
_ S 


« B 1 +6 2 6f 2 +C D a 2 «5 + +03(a'-a^ T ) 

w i 


> a 


i + > 
i + 


(2.3.26) 

In Figure 2.8 drag coefficient is shown as a function 
of angle of attack, o, for zero and maximum flap deflections. 
Lift coefficient is shown as a function of drag coefficient 
in Figure 2.9. 
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Figure 2.9 Lift Coefficient vs. Drag Coefficient 
on Wing 
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2.3.2 Tail 

During most of the flight regime under consideration 
the rotor downwash does not have a significant effect on 
flow over the tail. In hover or when the forward speed is 
low, the induced velocity is high. Under such conditions 
the downwash is pointing almost vertically downward and does 
not pass over the tail. When the forward speed is high and 
the transition is over the induced velocity is so low that 
it does not produce any significant effect on the tail. 

We use drag and lift values on the tail for zero ele- 
vator deflection. The direct effect of elevator deflection 
on vehicle lift and drag is neglected. The indirect ef- 
fect through change in angle of attack on the wing is im- 
portant. 

The angle of attack on the tail is 

aj=e-y+o^j (2.3.26) 


where c^.j is the angle of incidence of the tall. The lift 
and drag coefficient on the tail are assumed to have ex- 
pressions similar to wing, i.e., 


C, ( #t ) = a T a 
l T 1 


IT 


a T I — a S T 


a T a sT[? 

" [7 -stD 

_ a T a sT + a _j} 

’ tt ■ *st3 


°T > °sT 


°T K -°sT 

(2.3.27) 
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C D t s C D t +C D t a T°T 
T 'O T 1 


°Tl ± °sT 


= C D t +C D t a T a sT +e 3 x ^ a T" a sT^ °T > a sT 
1 0 1 1 1 


8 C D, +C D t a T a sT + e 3 T ^" a sT‘ a T^ a T < ‘“sT 
‘O T 1 1 


(2.3.28) 


The drag and the lift on the tail are given by 


4*1' v2s T c L T 

d t * r » v2s t c d t 


(2.3.29) 

(2.3.30) 


2.3.3 Fuselage 

The lift and drag on the fuselage depends on the shape 
of the fuselage in addition to the fuselage angle of attack. 
The fuselage is approximated by a long cylinder. The drag 
coefficient and lift coefficient can be written as ( HO- 1 ) 

C D = C D si n 3 1 ctp | +C Q (2.3.31) 

F F 1 F 0 

p 

C, = C, sin or cos or sgn (a p ) (2.3.32) 

l F l F 1 t ■ * r 

The lift coefficient and drag coefficient are shown in 
Figure 2.10 as a function of the fuselage angle of attack, 
ap , given by 

Op = 6-y (2.3.33) 


28 



Figure 2.10 Lift and Drag Coefficients 
on the Fuselage 

If Sp is the fuselage platform area, the lift and the 
drag can be written as, 

D F = \ pV 2 SpC Di _(op) (2.3.34) 

L F = \ pV^ SpC^^(ap) (2.3.35) 

2.3.4 Nacelles 

The nacelles are completely immersed in rotor downwash 
under all flight conditions. They are also approximated by 
long cylinders. The lift and drag coefficients of equations 
(2.3.31) and (2.3.32) can then be used. If is the 

platform area of the two nacelles, the drag and lift are 


29 



given by 


°NAC “ 7 p V R S NAC C D NAC ^“NAC^ (2.3.36) 

L NAC ' ? P V R S NAC C L NAC (a NAC^ (2.3.37) 


where V R is given by equation (2.3.4) and « NAC is the 
nacelle angle of attack. 

a NAC = n+Op-Ao (2.3.38) 
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CHAPTER III 

Tilt-Rotor VTOL Aircraft Guidance 
3.1 Introduction 

The requirement for a good guidance scheme for near 
terminal operations of a VTOL aircraft was discussed in 
Chapter I. The scheme should be as simple as possible so 
that a small airborne computer can handle it. We have 
adopted the perturbation guidance scheme because of the 
simple linear feedback law that it gives. This method 
consists of linearizing the nonlinear state equations about 
a nominal path. The feedback control gains are obtained 
by using this linearized system of equations and suitably 
chosen weighting matrices in the quadratic synthesis tech- 
nique . 

The equations of motion in the vertical plane are 
derived in Chapter II. They are, 


= F 1 ( V , y ; 0 , n » T , 6 f ) (3.1.1) 

V ^ = F 2 ( V ,y; 6 ,n »T,6f ) (3.1.2) 

= V sin y (3.1.3) 

d x 

= v cos Y (3.1.4) 
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The symbols have the same meaning as in Chapter II. 


Under almost all flying conditions it is possible to 
generate large control moments. In the helicopter mode 
pitching moment is produced by dual longitudinal cyclic 
pitch change, yawing moment by differential shaft tilt of 
the two rotors (or by differential longitudinal cyclic pitch 
change), and rolling moment by differential collective 
pitch change. In the cruise mode control moments can 
be produced using the regular aerodynamic surfaces or cyclic 
and collective pitch changes. Because of large available 
pitch moment generating capability it is possible to put 
a tight control on the pitch angle, under these conditions, 
the pitching motions occur in a frequency range very much 
higher than the frequency range of the translational motions 
governed by equations (3. 1.1-4). Therefore the transla- 
tional motion equations can be decoupled from the angular 
motion equations and 9 can be taken as a control variable 
for the translational motion. In this system, V, y, x 
and h are state variables and e, ii, T and 6f are control 
variables. Since none of the equations depend explicitly on 
t, we can use x as the independent variable. The equations 
of motion then become, 

dV _ F 1 

<J7 “ V COS y 

dx = F 2 

V Z COS y 


f 1 (3.1.5) 

f 2 (3.1.6) 


|jj- = tan y (3.1.7) 

*Later we shall use h as the independent variable during the 
vertical and near-verti cal part of the take-off. 
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If time, t is required it can be obtained by solving 


dt = 1 

cTx ” V cos y 


(3.1.8) 


3.2 Linearizing of the Equations of Motion 

The guidance scheme must work in the presence of steady 
or fluctuating winds. Therefore, it is necessary to distin- 
guish between air speed and ground speed (see Figure 3.1).* 
The air values are represented by subscript "a" and ground 
values by subscript "g". Nominal values are represented by 
a bar. In the presence of wind the lift and the drag do not 
act perpendicular and parallel to the ground velocity. The 
equations of motion must be modified to 


w HORIZONTAL WIND 



Figure 3.1 Ground and Air Speed in the 
Presence of Wind 


*Near the ground vertical winds are negligible compared to 
horizontal winds 
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where 


dx 2, ' ^ F i( v a * Y a* 0 » n,T,6f) cos ^ Y a" Y g ) } 

_F 2^a ,Y a * e * n ,T * sin (Y a "Yg)| / (VgCos Yg ) 


F 1 a cos ^Ya"Yg)~ F 2a sin (Ya"Y g )j / ( V g cos Y g ) 

(a. 2.1) 


dx IT? 


V COS y 

g g 


{ F 2a cos (Y a‘ Y g )+F la s1n(Y a' Y g>} 


dh _ . 

7x ■ tan Y g 


V a - V g ^ 2 + 2V gW COS Yg 


( 3 . 2 . 2 ) 

( 3 . 2 . 3 ) 


sin y 


= ll 

a V. 


sin Y , 


If equations governing nominal values of ground speed, 
flight path angle and vertical distance (where wind is zero) 
are subtracted from the above equations, we have 


d V dV 

3T - 3T ■ CF la cos (va-’ f g ) - F 2a s1n(l 'a- 1 'g )]/(V g cos V 

“ F 1 ( V g | Yg » 0 > T) » T 9 6 f ) / V g CO S Yg] 

" C F ia cos(l, a- lf g ) " F 2a sl ' n(Y a' Y 9 ):l/(v g cos Y g> 

F ig + F m In 

V os Y g v g cos Y g v g c ° 5 % 


34 


. ( F 1a cos (T a -Tg>- F 2 a s1 " 


V COS y 

g g 


+ ^la 


111 


v g cos Y g v g cos 9 g 


(3.2.4) 


The first term represents the effect of wind and the 
second term the deviation of the state and control variables 
from nominal values. Now 


IlL 


' lg 

V g C0S Y g V cos y 


- f -f - 

" f ig f ig 


„ 9 f 


3V. 


l! 


3f ln 

6 V + -r— i-3- 

9 3y g 


af 


6 Y + 

y g 36 

g 


la. 


66 




9n 


3f 


6ti + 


IT 


Is. 


3f 


6T + 


3(6 


& 


6(6f) , 


(3.2.5) 


for small changes. 


The first term in (3.2.4) can be simplified as, 


F la cos (Y a -y n )-F ?a sin ( Ya - Yn )-F 


a r g / '2a 


a T g' 'lg 


V cos y. 

g g 


F 1 a~ F lg~ F 2a (Y a~ Y q ) 
V„ cos Y , 


(3.2.6) 


g 

1 


V_COS Y- 

g g 


g 

3F lg 8V a 3F lg 3y a 

3V _ 3w 3y 3w 

g g 


9 r 

2g 3w 


,1 

“-I w 
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where w is the wind velocity. Since 


F ig ■ 

f ig v g 

3F lg 

w~ 

. !> 
3V„ 

g 

g 

Ilia 

. ilia 

8Y 9 

3y g 


ig 


ig g 


(3. 


(3. 


Therefore , 


F la c° s (T r,-Y q )-F 2a s1n (Y a -T q )-F 1q 


V g cos y g 


+ (^ • f ig tan T g' Vs) $ 


w (3. 


hence equation (3.2.4) becomes. 


d 6V 


ax 


3 f 

IV 


3 f ■ 


3 f 




+ ~- 9 - 6n + 6T + <5 ( 6 f ) 


3n 


3T 




11 

aw 


+ (p4. 

\3y„ 


f lg tan Y g" V g f 2g 


) ^} W 
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2.7) 
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2.9) 
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3 f p n 

^ sv 9 ^ 


8 f p n 8 f p 

® Y g + Te -1 se + ITT 3, 


3f 2a 9f ?n f/ 9f ?n 2f, ni aV a 

+ vr* 4T + ist 2 5({f > + ((■ sv^ + -ir 1 )^- 


♦ (1-21 


f 2g tan Y g + r^jlw 1 } " (3 ‘ 2 • 1 1 ) 


* SeC Y g % 


(3.2.12, 


All the coefficients refer to nominal path. 

It is necessary to evaluate 9V a /9x and g Y /9x as a 
function of V g and Yg (see Figure 3.1). 

V a = (Vg+w 2 +2V g w cos Y g ) 1/2 (3.2.13) 

V 1 +& cos y + ^ -V 

3V V -VI 9 V g 9 9 

a. = a q _ *- a v q J 


V 

= ^ ?“ cos Yq = COS y (3.2.14) 


.-1 V in 


Y a = tan ' -2— 

a V cos y +w 

9 9 


(3.2.15) 


v „ s in Y, 


!li . Llm (V 3 C0S V w) 

9w " w+0 77. : rj” 

w+0 T (V q S1n Y o ) 

(V g cos Yg+w ) 2 
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V s i n y 

— *_ 2 _ . — La_ — „ — 

V cos y n + V 2 sin y n 

g g g g 


sin Y ( 

rr" 


(3.2.16) 


The subscript g is removed now with the understanding 
that all values refer to nominal ground values. Arranging it 


in the state space form, 

{jf = F ( x ) Z + 6 ( x ) U 


where 
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9 V 
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" 9f l 
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r (x)w (3.2.17) 
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g 2e 

g 2n 

9 2T 

9 26f 
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U 


6n 

6T 

6 ( «f L 


r(x) = 


(f 11 +f 1 /V)cos Y- ^ (f "i 2 "^i tan Y -f 2 v ) 

• / f -i 

( f 2i +2 f 2 /V ) c° s y- [f 22 -f 2 tar\ y+ y- 


J 

(3.2.18) 


There are four control variables: Incremental pitch 
angle, 66 (changed by elevator deflection or rotor cyclic 
pitch control), incremental rotor shaft tilt, 6n, incre- 
mental thrust, 6T (varied using collective pitch control), 
and incremental flap deflection, 6(6f). 

3.3 Quadratic Synthesis 

For the system given by equation (3.2.17) we use the 
quadratic synthesis method to minimize 
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J = | Z T (x f )S Xf Z(x f ) 



0 


( 3 . 3 . 1 ) 

where A and B may be constant or functions of the independent 
variable. One possible choice of S x ^» A and B is suggested 
by Bryson and Ho (BR-1). This choice uses a constant A 
and B and is good as a starting value in most cases. However, 
in this system this may not be true because of the change 
in configuration from helicopter to cruise mode and vice 
versa. In other words the system definition matrices, 

F, G and r vary with x considerably. The system is very 
sensitive to variations in n when the speed is low. If 
constant values are taken from matrices A and B, the gains 
tend to be undesirably high at low speeds. 

While it is clear that it is desirable to use x-varying 
matrices A and B, the choice of x-dependence is not clear. 

One reasonable way, we think, would be to choose 

A~^ (x) ^ max Z(x)Z T (x) 


For this case these equations were found to be very 
useful in making an initial guess. To simplify things, we 
used B depending linearly on V. 
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The solution to the optimization problem (equations 
(3.2.17) and (3.3.1)) requires the solution to the well- 
known Riccati equation 

|| = -SF-F T S-A+SGB" 1 G T S S(x f ) = $ x ^ 

(3.3.3) 

The equation can be integrated backward starting from 
the known final condition. The control can, then, be found 
as a function of the state variables, 

U = -B -1 G T SZ = CZ (3.3.4) 

3.4 Control in the Presence of Wind 

In the presence of wind, it is helpful to feed back 
an estimated value of the wind with a gain. To derive an 
equation for this gain, a model for the wind velocity is 
needed. A reasonable model is 

w = Ew + white noise (3.4.1) 


Adjoining this equation to the other state equations 
(3.2.17) we have: 


r*i 



f---l ♦ 

U-l 


O 

m 

L w J L 0 J 


or, 
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(3.4.2) 


It can be shown (BR-2) that the gains on 6V, 6y» <5y a' 1 ' 6 
the same as when w = 0, and that the gains on the estimated 
wind w are obtained by solving the following differential 
equation for R 


R = -Sr-RE-F T R+SGB _1 G T R 


R(x f ) = 0 

(3.4.3) 


where , 



and C w = -B-VR (3.4.4) 


3.5 Estimation of the Wind Velocity 

To control the aircraft effectively in the presence of 
wind, it is necessary to obtain an estimate of the wind 
velocity. It could be found by direct measurement or from 
some other external source. However, if this is not possible, 
the wind velocity may be estimated using measurements of 
aircraft motion. This method of estimation is described 
here. 


The state equations and the equation governing the 
average wind velocity are 
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Z ■ F ( x ) Z + G(x)U + r(x)w 


(3.2.17) 


w = Ew (3.4.1) 

Take Z as the measurement (assuming Z is available). 
Then the estimated value, w, of w can be obtained by 
• A 

w = Ew + K(x) [Z-Z] (3.5.1) 

From equation (3.2.17), 

/V 

Z = FZ + GU + rw (3.5.2) 

Substituting equations (3.2.17) and (3.5.2) in equation 
(3.5.1) we have, 

■ 

w = E w + K(x)r(x)[w-w] (3.5.3) 

Subtracting equation (3.4.1) from equation (3.5.3) 


^ (w-w) = E(w-w) - K(x)r(x) (w-w) (3.5.4) 

Defining the error e as, 

e = w-w (3.5.5) 

gf = ^E-K(x) r ( x)j e (3.5.6) 

By a suitable choice of the x-varying gain K(x) the 
characteristic root of (3.4.6) can be placed anywhere one 
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TWO IMPLEMENTATIONS ARE POSSIBLE 


(i) if Z is available. 



(ii) if Z is not available 



K(x) 



Figure 3.2 Estimation of Wind Velocity 
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desires, for all x. 


Equation (3.5.3) can be put in another form using only Z 
and not Z : . 

w-KZ-KZ = Ew+Krw-Krw-KZ-K[FZ+GU+rw] 


i . e . . 


defining 


fx <"- KZ) 


y = w-KZ, 


= (E-Kr) 


jw-KzJ -KGU-|KR+ME-Kr)l<Jz 

(3.5.7) 


4 — y = (E-Kr)y - KGU - ( Kr+fc- ( E-Kr) K) Z 

(3.5.8) 

v This is another way the observer can be implemented. The 
gain K(x) can be chosen as before. The block diagram for 
two possible observers is shown in Figure 3.2. 


3.6 RMS Response in Presence of Fluctuating Mind 

A tilt-rotor aircraft may encounter fluctuating wind 
velocity along its flight path because of variation in 
wind velocity with (a) distance, and (b) time. When 
the correlation distance is small, it is not possible to 
estimate the wind velocity accurately. If the correlation 
distance is infinitesimally small, i.e., the wind is "white", 
the gains on the estimated wind do not change the response 
since it is impossible to predict the future behavior of 
the wind from past information. Substituting for U from 
equation (3.3.4) in equation (3.2.17), we have 
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-j-j- = (F+GC)Z+rw . (3.6.1) 

If w is assumed white, with time correlation 

E[w(t)w(x)] = Q ( t ) 6 ( t-x ) (3.6.2) 

2(x) = E(Z (x)Z T (x) ) (3.6.3) 

it can be shown that (BR-1), 

= ( F+GC ) Z + Z ( F+6C ) T + rQ(x)r T (3.6.4) 

where , 

E(w(x)w(x ' ) ) = Q(x) 6 (x-x* ) (3.6.5) 


Since, 

dt = 1 

H3T - V cos y 

equation (3.6.1) can also be written as 

^ = ( F+GC ) Z V cos y + rw V cos Y (3.6.6) 

Hence 

= ( F+GC ) V cos yl+l V cos Y (F+GC) T +rQ(t)r T V 2 cos 2 Y 

(3.6. 7) 

Changing back to x as the independent variable, 

4| = (F+GC)Z+Z(F+GC) T +r Q(t)V cos y r T 



Comparing (3.6.4) and (3.6.8), we have 


Q ( x ) = Q ( t ) V cos y (3.6.9) 

Equation (3.6.9) is a general equation for changing 
the independent variable. Equation (3.6.4) can be solved 
using equation (3.6.9) and the initial condition, 

2(x 0 ) = Ejz(x Q )Z T (x 0 )j (3.6.10) 

3.7 Summary 

This chapter presented a general perturbation guidance 
technique for a tilting proprotor aircraft. This scheme 
can be simplified in most cases. For instance, it may not 
be necessary to use all four control variables or it may be 
unnecessary to feed back all state variables to all controls. 
The next two chapters consider take-off and landing tra- 
jectories for the Bell Model 266 and show how simp! i cati ons 
can be made in those cases. 
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CHAPTER IV 

Longitudinal Guidance During Take-off 
4.1 Introduction 

The general method for VTOL aircraft guidance, pre- 
sented in the last chapter is modified and simplified to 
develop a guidance scheme for the take-off phase of 
the flight. The take-off involves transition from a low- 
speed helicopter mode to a high-speed cruise mode, in which 
the rotors are used as propellers. The techniques applied 
here can be used for any take-off trajectory. However, 
a particular flight path has been chosen for illustration. 

The selected flight path is shown in Figures 4.1a 
and 4.1b. The aircraft starts vertically from zero speed. 



Figure 4.1b Terminology Used for Take-Off 
Trajectory 

The flight path angle is maintained at u/2 until a speed of 
2.85 m sec"^ is reached. Then the aircraft starts the 
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Function of x for Nominal Trajectory 





transition at a constant normal acceleration of -1.6 m sec 
The thrust is fixed at its maximum available value at all 
times. The flaps and aileron-flaps are deflected as much 
as possible (57 degrees) before starting, to reduce 
download on the wings, and are retracted as the aircraft 
gains speed. The pitch angle is chosen so that the wings 
are loaded as much as possible with a 20% margin for stall. 

We specify maximum and minimum allowable values of the 
pitch angle. Figure 4.2 shows nominal values of flight 
path angle, velocity, pitch angle and rotor shaft deflec- 
tion angle. 

4.2 Simplifying the Model 

The problem of guidance during take-off has control 
variable constraints. The thrust is fixed at the maximum 
value; hence positive values of 6T are not feasible. 

There are similar constraints on flap deflection and pitch 
angle during a portion of the take-off trajectory. One 
way to use quadratic synthesis is to use less (more) than the 
maximum (minimum) values of these control variables during 
the nominal flight leaving a certain margin for feedback 
control. However, in this case, there are four control 
variables and it is not necessary to use all of them to 
achieve the desired goal. We decided to use thrust and 
flap deflection as open-loop controls and rotor shaft tilt 
and pitch angles as feedback controls. The installed 
power and rotor size in a tilt- rotor VTOL aircraft are 
dictated by the take-off requirements; hence using any- 
thing but the maximum thrust during take-off is economically 
undesirable. The flap deflection is not a very effective 
control variable. The linearized equations in the pre- 
sence of wind (equations (3.4.2)) simplify to 
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(4.2.1) 


IT = F w (x)Z w + G T0 (x) U T0 

where F , Z , r and w are defined in equations (3.2.18) and 
w w 

(3.4.2) and 6 jq and Ujq are 


G 


TO 


g ln 9 le 

g 2n g 2e 
0 0 

0 0 


u 


TO 


Sn 

60 

w 


(4.2.2) 


The drag on the rotors and the nacelles is small along 
most of the flight path and is dropped. 


4.3 Choice of Weighting Matrices 

Tentative specifications for matrices A, B, and S x 
in equation (3.3.1) were made using equations (3.3.2) and 
(3.3.3) as guidelines. A simplified linear V-dependence is 
used. Final values of these matrices were obtained by 
cut and try. They are* 


.00021 


A = 


0 

0 

0 


0 0 0 

.04 0 0 

0 .00044 cos Z y 0 

0 0 0 


♦Units in meters, sec, rad. 
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.14 V 




,-l 


3^ V>24 . 5 m sec" 1 
0 J V <24 . 5 m sec" 1 


S 


f 


.069 1.2 .02 

1.2 170.0 2.9 

.02 2.9 .063 

-.0026 1.1 .013 


-.0026 
1.1 
.013 
.048 _ 


(4.3.1 ) 

2 

The Agj term in matrix A varies as cos y because we 
are interested in minimizing the square of the distance 
normal to the flight path while the state variable Zw 3 is 
the vertical distance. It is clear from Section 3.4 and 
the fact that the wind velocity w, is uncontrollable, that 
the weight A^ 4 makes no difference. At low speeds, the 
angle of attack of the wing, and hence the pitch angle, does 
not affect the forces very much. Below a speed of 24.5 
m sec" 1 , we put (B‘ 1 ) 22 equal to zero to reflect this. The 
matrix is chosen as the steady state Riccati gain matrix 
for a constant coefficient system with state and weighting 
matices corresponding to the terminal time.* This ensures 
that no gain is abnormal near the terminal condition. In 
cases where the process continues beyond the terminal design 
time, we consider this a useful algorithm for picking the 

Riccati gain matrix at the terminal point. 


* 

In this case S v is computed using w = - . 0001w+"wh i te"noi se 
x f 

to avoid the problem of zero eigenvalue in the eigenvector 
decomposition technique. 
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4.4 Control Gains and Aircraft Performance in the Presence 

of Wind 

The gains on rotor shaft tilt and pitch angle are shown 
as a function of the horizontal distance x in Figures 4.3 and 
4.4. These gains are for aircraft speeds greater than 2.85 
m sec"^. During the vertical flight immediately after take- 
off, the horizontal distance does not change and it cannot 
be used as the independent variable. This problem is dis- 
cussed later in Section 4.5. 

All gains are negative or zero throughout the flight 
path. The gains on pitch angle are zero for x < 64 m 
since we put ( B” ^ ) 2 2 e 9 ua ^ t0 zero for x this ran 9 e * 

Notice that most of these gain programs can be approximated 
by combinations of constant values and ramp functions. 

Controlled tilt-rotor aircraft flight paths in the 
presence of wind or initial disturbances were computed 
using the nonlinear equations. At any point along the 
flight path, incremental pitch angle and rotor shaft tilt 
are computed using deviations in aircraft velocity, flight 
path angle, vertical distance from the nominal trajectory, 
and estimated wind velocity. A continuous implementation 
is assumed. 

Figure 4.5 shows the flight path for the tilt-rotor 
VTOL aircraft with initial errors in flight path angle 
(-5.7°) and distance away from the nominal trajectory 
(vertical distance = -5m). The errors die out quite quickly. 
The normal distance of this trajectory from the nominal one 
is almost zero at the final horizontal position. With no 
control there will be a terminal error of 10m. The 
figure also shows variations of 6e and 6n. By the time the 
pitch angle control comes into action the errors are already 
small. The maximum value of 6.n is less than 6° and <se is 
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Figure 4.4 Gains on Pitch Angle 
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Figure 4.6 Trajectory and Control Effort in the Presence of 5 m sec" 
Head Wind [Wind Velocity Estimated]. 



always small. 

A trajectory in the presence of a steady head wind 
of 5.0 m sec' 1 is shown in Figure 4.6. Initially this 
trajectory falls below the nominal trajectory. The 
error in position starts from a low value, reaches a maximum 
value of 1.4 m at horizontal distance 300 m and then de- 
creases again. The final error is less than .4 m. In the 
absence of feedback control it would be 35 m. The devia- 
tions in rotor shaft tilt angle and pitch angle reach their 
maximum values at 200 m. By this point the wind velocity 
is estimated within 5% and there is enough 66 and 6n to 
counterbalance the effect of wind as much as possible. This 
required amount of 66 and 6n decreases with increasing speed. 

Using control gains computed above many other flight 
paths were computed. The flight paths obtained by assuming 
perfect knowledge of wind were superior to the ones with a 
wind estimator. Feeding back estimated or otherwise measured 
wind velocity provides a considerable improvement. 

Using the linearized system of equations, the root mean 
square (RMS) response is computed in the presence of random 
wind. Figures 4.7 and 4.8 show RMS values of state and con- 
trol variables in the presence of a random wind with RMS 
value of 3.0 m sec' 1 and correlation time of .30 sec 
(modelled "white"). All state variables are assumed to be 
known exactly. The RMS values of initial errors in velo- 
city, flight path angle and height were taken as 1.0 m sec' 1 , 
.01 rad, and 3.0 m respectively. 

Figures 4.9 and 4.10 show RMS values of state and 
control variables in the presence of wind consisting of two 
components, (a) a random wind as in the last case, and (b) 
a steady wind with RMS value 3.0 m sec" 1 . The initial RMS 
errors are the same as in the last case. Initially the 
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Figure 4.7 RMS State In the Presence of Random Wind ["white" with Power 
Spectral Density 5.4 m 2 sec _1 ] 
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Fiqure 4.8 RMS Control in the Presence of Random Wind 

2 

["white" with Power Spectral Density 5.4 m sec 
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Figure 4.10 RMS Control in the Presence of Random Wind of Figure 4.9 




errors increase and then they decrease. Within 100 m from 
the starting position the RMS errors reach almost a con- 
stant value. The RMS error in velocity is about 1.0 m sec”^ 
and in height is about 1.0 m. These errors are more than 
the corresponding errors in the last case where there is 
no steady wind. The RMS 6n, the change in rotor shaft tilt 
angle, is always less than 6° and the RMS 6e never exceeds 
1.5°. Less control is required in the absence of this 
steady wind. The errors in state variables and the required 
amount of control look quite reasonable. 

4.5 Guidance at Low Speed 

Immediately at take-off the flight path is vertical and 
horizontal distance is not a useful independent variable. 
Furthermore the assumption concerning the decoupling of the 
longitudinal and lateral vehicle motions is not very good. 

One way to deal with these complications is to use a 
modified hover autopilot until the aircraft reaches a cer- 
tain speed. Since commercial VTOL aircraft will probably 
have an autopilot hover mode, it should be easy to modify 
this mode so that it can be used at low speeds also. Hover 
autopilots for VTOL aircraft have been designed for example 
by Bryson and Hall (BR-2). 

Even if the assumption of decoupling is considered satis- 
factory for this purpose there is another problem. At low 
speeds even when the flight path is not exactly vertical 
the gains and nominal values of state and control variables 
change rapidly as a function of the horizontal distance, x. 
Errors in estimation of x could cause some trouble. Under 
these circumstances it is better to use the height h as the 
independent variable during the initial part of the trajectory. 


64 


- 2.0 



65 


Figure 4.11 Gains on Rotor Shaft Tilt with Height as the Independent Variable 







To do this it is necessary to relinearize the equations of 
motion. 

Altitude,- h, is used as the independent variable until 
the flight path angle is 45° and thereafter it is changed 
to x. The flight path angle is 45° at horizontal distance 
27 m and height 41.7 m. The gains are shown in Figure 4.11 
for 0 < h < 41.7 m. In these computations the weighting 
matrices were not modified to account for a change in the 
independent variable. The gains do not change as rapidly as 
they did before. 


4.6 Summary 

Table IV-1 shows the control system performance with 
initial disturbances. The initial disturbance is -5.0 m in 
height and -.1 rad in flight path angle. The control system 
cuts down the error by a factor of about 100 while using a 
nominal amount of control. 

Table I V - 2 compares the control system under different 
conditions for a 5.0 m sec~^ steady head wind. In the ab- 
sence of control, the average error (deviation normal to 
the nominal trajectory) is 24 m and the error at the be- 
ginning of cruise is 35 m. If the wind velocity is not 
estimated or measured the average error is 5.0 m and the 
error at the beginning of cruise is 2.6 m with the control 
system. The change in rotor shaft tilt angle averages 
6.7° and reaches a maximum value of 11°. If the wind 
velocity is estimated the average error is 0.8 m and the error 
at the beginning of cruise is .4 m. The change in rotor 
shaft tilt angle still averages 6.7° and has a maximum of 
11°. The performance is best when wind velocity is exactly 
known followed by the case where it is estimated. 
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Disturbances [e. =-5.0 m, e -.1 rad] 










<D 




"O 


c 


03 


cr 


<o 


M— 


O 


03 


to 3 


0) i— 


3 03 


r— > 


03 


> E 


3 


E E 


3 


E X 


*1— 03 


X E 


03 


E <D 


-C 


03 +-> 


.C 


+-> to 


•i— 

X 

C 


o to 

03 

■r— 

O 

-o -C 

c 

0) 1 — 

03 

to 

4-> 

o 

to 

Q- • 

•r— 

E c 

XJ 

•1“ o 


•1— 

r~ 

4-> +-> 

03 

C 03 

■M 

•i- 3 

C 

03 cr 

o 

U 03 

N 


•t— 

to S- 


C 0 3 

o 

O 03 

-C 

a c 


•r— 


03 « — 

03 

a 

> 

to o 

o 

03 C 


* 

03 

03 

03 

03 x: 

03 

-M 


03 

03 

-C C 

> 

h“ *r- 

<£ 

■K 

4- 


68 


Table I V - 2 Comparison of Different Control Laws in the 
Presence of 5.0 m sec"^ Head Wind 











The x-varying control gains were approximated using a 
combination of constant values and ramp functions. The 
approximated gains are shown with dotted lines in Figures 
4.3 and 4.4. This approximation reduces the storage re- 
quirement and computation time for an airborne computer 
considerably. Referring to Tables I V- 1 and IV-2, we see 
that the approximate gains deteriorate the performance 
slightly. The approximated gains should be used in 
actual implementation. 
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CHAPTER V 

Longitudinal Guidance During Landing 

5.1 Introduction 

This chapter deals with the control gains and the 
guidance scheme for a tilting proprotor VTOL aircraft during 
approach and landing. The technique is similar to the take- 
off case. 

The nominal trajectory considered here is a straight 
line with a flight path angle of -5.7°. The initial 
aircraft velocity is 66.8 m sec" 1 and the initial flap 
deflection is 14 degrees. To start with the rotors are 
unloaded and rotated at a rate of 5.7° sec' 1 until the 
rotor shaft angle reaches 101 degrees. The flaps and the 
aileron-flaps are also turned at 5.7° sec" 1 rate until 
they are deflected 57 degrees. The rotors start producing 
thrust when the rotor shaft angle exceeds 90° so that a 
component of the thrust decelerates the aircraft. The pitch 
angle is fixed at -6° during powered flight.* The 
aircraft trajectory is shown in Figure 5.1. Nominal values 
of velocity, thrust, pitch angle and rotor shaft deflec- 
tion angle are shown in Figure 5.2. 

5.2 Simplifying the Model 

The linearized equations of Chapter III have been 
simplified as in the case of take-off guidance with one 
exception. During take-off the rotor thrust was set at its 
maximum available value at all times. However, during 
landing the thrust is less than the maximum value and can 
be used as a control variable. We use rotor shaft tilt 
and rotor thrust as two control variables. The simplified 
and linearized equations are 


*0therwi se the pitch angle is determined by the requirement 
of zero normal acceleration. 
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Figure 5.2 Velocity, Rotor Shaft Tilt and Rotor Thrust 


for the Nominal Trajectory 
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(5.2.1 ) 


dZ 


IT 


-- F w (x)Z W + G L<*> U L 


where 


and 


r 


G L = { 


g 1n 

g 2n 

0 


g lT 

g 2T 

0 




J 



<$n 

6T 


(5.2.2) 


5.3 Choice of Weighting Matrices 

The technique for choosing weighting matrices is very 

similar to the take-off case and is described in Section 

4.3. The final values of A, B and S v are (units meters, 

x f 

sec, rad) 


.0011 0 


0 0 


A = 


0 

0 


.04 0 0 
0 .00044 cos^y 0 


L 0 0 


0 0 J 


".082V 0 

_ 0 2x10 7 V- 


and , 
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-.0036 

.0044 

.0072 
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i 

-.06 

-.057 

240 J 



(5.3.1) 

5.4 Gains and Flight Paths in the Presence of Wind 

Figures 5.3 and 5.4 show gains on rotor shaft tilt and 
thrust as a function of the horizontal distance. Near the 
terminal point, the optimal gains had wide fluctuations 
in spite of choosing S x ^ given above. The gains have been 
smoothed in this region. 

Notice that gains on rotor shaft tilt are zero for 
x < 930 m. For this range of x the nominal thrust is zero. 
Hence changing the rotor shaft tilt produced no change in 
forces . 

Unlike the take-off guidance some gains in this case are 
positive. There are two main reasons for this: (a) Since 
the landing trajectory is a straight line the velocity- 
dependent centrifugal force is absent, (b) during a major 
portion of the flight the rotor shaft is tilted more than 
90°; under this condition increasing the rotor shaft tilt 
produces many effects which are similar to the effects pro- 
duced by decreasing the rotor shaft tilt if it were less than 
90°. 


In the absence of the centrifugal force, it is much 
easier to explain the nature of the control gains. If the 
rotor shaft tilt exceeds 90°, increasing it decreases the 
force along the velocity vector and perpendicular to the 

*In this case S x ^. is computed using w s - . 0001w+"whi te" noise 
to avoid the problem of zero eigenvalue in the eigenvector 
decomposition technique. 
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Figure 5.3 Gains on Rotor Thrust 
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Figure 5.4 Gains on Rotor Shaft Tilt Angle 
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velocity vector. This tends to reduce velocity, flight path 
angle and height. Hence C^y, C and are positive or 
zero . 

The gains on rotor thrust are never zero. It is 
assumed that it is possible to produce small negative 
thrusts. During the initial part of the trajectory 
when the rotor tilt is less than 90°, an increase in thrust 
increases the force along the velocity and perpendicular 
to the velocity. Thus Cyy, and C-^ are negative. As 
the rotor tilt exceeds 90°, Cyy tends to change sign. 

A positive wind velocity (head-wind) has two effects. 

It decreases the angle of attack (because the nominal flight 
path angle is negative) and increases the air speed. It is 
not straight forward to predict the signs of C^ w and C^ w . 

Figure 5.5 shows the errors in velocity and the distance 
perpendicular to the nominal flight path for an initial error 
of -0.1 radian in flight path angle and -5 m in vertical 
distance. This is a large error since it increases the 
angle of attack by 5.7°. The error in velocity increases 
to about 2 m sec"^ and then decreases. The distance normal 
to the nominal flight path increases from -4.9 m to -12 m 
The aircraft reaches the terminal point with a small error 
in position and a small velocity in the forward direction. 

The figure also shows 6r> and 6T. The change in thrust 
starts from a high positive value and soon decreases to a 
low value. The <5n never exceeds 2.5°. With no control the 
velocity would be 5.1 m sec"^ and distance normal to 
nominal trajectory would be -9.5 m at the terminal point. 
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Figure 5.5 Errors and Control Effort with Initial Errors 
In Position and Flight Path Angle 


78 


2 


Error in 
Velocity, 
[m sec-'] 





Figure 5.6 Errors and Control Effort in the Presence 
of Head Wind 
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Figure 5.6 is a similar diagram when the aircraft 
flies in a constant head-wind of 5.0 m sec" 1 with no initial 
errors. In the absence of feedback control the aircraft 
never reaches the terminal point. The forward velocity is 
zero 165 m away from this point and the position error 
(distance normal to the nominal trajectory) is -80 m. With 
feedback control the error in velocity never exceeds 
6 m sec" 1 and distance normal to the nominal trajectory is 
always less than 6.0 m. The aircraft comes to rest 3 m away 
from the terminal point and the error in position is 2.8 m. 

In this computation all state variables are assumed to be 
known exactly except the wind velocity. The wind velocity is 
estimated using the technique of Section 3.5 with a filter 
of time constant 4 sec. The maximum value of |6n! is 5.7° 
and the maximum value of 1 6 T | is about 2000 N. As in the case 
of take-off the controlled trajectories were computed 
using the nonlinear equations. 

Figures 5.7 and 5.8 show RMS values of state and control 
variables as a function of the horizontal distance, x, in the 
presence of a horizontal wind with two components: (a) a 
steady wind with 3.0 m sec" 1 RMS value plus (b) a random 
wind with 3.0 m sec” 1 RMS value and .30 sec correlation time. 
RMS values of initial errors in velocity, flight path angle 
and vertical distance are taken as 1.0 m sec" 1 , .01 rad and 
3.0 m respectively. The linearized equations are used for 
this computation since the use of nonlinear equations makes it 
necessary to use the expensive Monte-Carlo type simulation. 

The RMS errors in velocity and position increase 
initially, decrease and then increase again near the terminal 
point. These errors are always reasonably small. The error 
in flight path angle stays small rising to a very high value 
when the aircraft speed is small. The control effort re- 
quired to obtain this performance is quite reasonable. 
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RMS State in the Presence of Random Wind, w 
[w=w-|+W 2 » w-j=random bias (RMS 3 msec"b* 

W 2 =white noise, power spectral density=5.4 m^sec"^] 
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5.5 Summary 

Most of the control gains obtained in this chapter can 
be approximated by a combination of constant values and ramp 
functions. It is necessary to do this to save computation 
time and required memory on an airborne computer. The ap- 
proximated gains are shown in dotted lines in Figures 5.3 
and 5.4. 


The performance of the optimal and the approximated 
gains in the presence of initial disturbances in velocity 
and flight path angle is compared in Table V-l. The tilting 
proprotor VTOL aircraft has about the same terminal and 
average error with optimal and with approximate gains. There 
is a considerable improvement especially in terminal velocity 
and position through using the feedback control. 

Table V-2 shows the VTOL behavior in the presence of 
5.0 m sec ^ head wind. There is a significant improvement 
in performance when there is a feedback on wind velocity. 

The use of approximate gains deteriorates the performance 
slightly. With no feedback control the errors are so large 
that it seems impossible to land the aircraft. 
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ompanson of Different Controls in the Presence of 5.0 m sec 




CHAPTER VI 
Lateral Guidance 


6.1 Equations of Motion 

The lateral motions (sideslip, e, roll angle, <j> , and 
yaw angle, i|>) of a tilt-rotor aircraft are nearly uncoupled 
from the longitudinal motions just as in the case of 
conventional aircraft. In general these three motions are 
coupled. In the absence of wind, these equations for small 
13, <f> and ip , are the same as for a conventional aircraft 
(BL-1 ) , 
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( 6 . 1 . 1 ) 


Def i ni ng , 





z 


L 




( 6 . 1 . 2 ) 


equations (6.1.1) can be written in another form in the 
presence of a lateral wind w (see Figure 6.1). 
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Figure 6.1 Lateral Motions of Tilt Rotor Aircraft 
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(6.1.3) 
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(6.1.4) 

Y is the control side force (if there is one) and N 
and L are control moments about the yaw and roll axes 
respectively. In a tilt-rotor aircraft direct side force 
could be produced by a simultaneous lateral cyclic pitch 
change in both rotors. However to avoid this additional 
complication in design, we shall instead use differential 
collective pitch to give a control rolling moment. 


The lateral distance away from the nominal flight path, 
y, is given by 


y = V e 


(6.1.5) 


Yaw moment is produced by differential rotor shaft tilt 
or differential cyclic pitch variation in the longitudinal 
direction. Large moments can be generated quickly because 
of high available thrust. Thus the roll and yaw angles can 
be changed much faster than the lateral deviation, y. In 
other words e,y motions occur at lower frequencies (or larger 
time constants) than the <f >,t|/ motions. Therefore, we can 
approximately decouple the lateral equations into two systems. 


.. . ? • 2/ W V\ 

^ + C 3 ^ + t*> o ^ + C 4 4 > = b 2 N+oj q ^ G+ i 7 JC *) 

w 

*c^ + C 5 4>-C 7 ip+C 6 Tp = b 3 L+C 6 (e+y^) 


(6.1.6) 


and. 
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(6.1.7) 


e + C i e = (C-j+C-j )\p+C24> +c i y 
y = V e 

In (6.1.7) we may treat \p and <p as control variables. 

The lateral wind velocity can be approximated by a first 
order system excited by a white noise, i.e., 

tw +w = n n(t) = 0 and n(t)n(t) = Q6(t-x) 

y 7 (6.1.8) 

6.2 Outer Loop Control Gains 

The outer loop of the lateral guidance scheme is de- 
fined by equations (6.1.7) and (6.1.8) in which a linear 
combination of ip and <j> , <j>+(C-j+C^ )/C 2 </> , is the control vari- 
able. For automatic take-offs and landings, it is necessary 
to control the lateral position. Writing these equations 
in terms of the lateral position away from the nominal 
flight path. 

y+C-jy = V(C 1 +C]) 4 ,+VC 2 <(.+C 1 w y (6.2.1) 

T "y +W y " " 

Let us choose 

( C +C 1 ) 

4> + — C" 2 ' * = -k-|y-l < 2-y“ l< 3 w y (6.2.2) 

Then we have 

y+C-,y = VC 2 (-k iy -k 2 y-k 3 W y )+C lWy (6.2.3) 


TW y+ w y = n 

The first of equations (6.2.3) is a second order system 
forced by w y . A good choice for control gain k 3 is obtained 
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by putting the coefficient of w y in this equation equal to 
zero, i .e. , 

C i c y g 

k 3 = it : = - rr (6.2.4) 

y <fr 

There are several methods for choosing the other gains. 
One method is to choose the gains in such a way that it be- 
comes a constant coefficient second order system. If w 
and are the desired natural frequency and damping ratio 
respecti vely , 



k 


2 



e 1 
v 

r 

♦ 


( 6 . 2 . 5 ) 


Expressions for C v , C w and C 

J 8 y <t y ip 

C is the coefficient governing the rate of change of 
side force with sideslip. It is given by (BL-1) 


C 


y 


e 




(6.2.6) 


is almost constant and is usually negative. Its 
value for the tilt rotor aircraft under consideration is not 
available. A value of -.6 was estimated. 


If e is the pitch angle the stability 
and C are, 


derivatives C 

y 4> 


C 


y 


<t> 


= cos e 

sq d 


( 6 . 2 . 7 ) 


\ s1 " 9 
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Let us pick control gains to obtain a natural frequency 
of 0.25 rad sec~^.and a damping ratio of 0.707 


k 


1 


(•25) 2 
Vg cos e 


V 



.0064 
cos e 


rad m 


k 


2 


= .3535 + ( " • 6 ) xSq d 1 

g cos 0 mg cos e V 


■t 


.036 
COS 0 


l.lxl0~ 4 V 
COS 0 


rad sec m 


-1 


k 


3 


1.1xl0~ 4 y 

cos 0 


rad sec m 


( 6 . 2 . 8 ) 


During a major portion of the take-off trajectory the 
pitch angle is 20° and for the landing trajectory is -5.7°. 
These angles are used in expressions (6.2.8). The con- 
troller gains for take-off and landing are shown in Table 
VI-1 . 


The small angle approximation for sideslip is not good 
at low speeds. This problem is reduced by writing the 
equations in the form (6.2.1) and (6.1.6) and redefining some 
stability derivatives. Though the sideslip angle may be 
large at low speeds, the lateral displacement and lateral 
velocity are usually small. However, the decoupling assump- 
tion is at best approximate when the speed is very low. 


6.3 Inner Loop Control Gains 

Equations (6.1.6) and (6.1.8) define the equations 
governing the inner loop dynamics for the lateral motion. 
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GAIN 


TAKE-OFF 


LANDING 


K 1 

.0066 

.0065 

rad m" 1 



k 2 

rad m"^sec 

.037 - 1.1 x 10" 4 V 

.036 - 1.1 x 10“ 4 V 

K 3 

rad nT^sec 

1 

O 

X 

• 

X 

o 

1 

-p* 

< 


Table VI-1 Outer Loop Control Gains for 
Take-off and Landing 


92 





These differential equations for <j> and t|» are coupled and 
are forced by e and w y . The coefficient C 4 is usually 
small and the coupling is approximately one way. 

The required value of $ + ( (C-j+C^ ) / C 2 ) ^ at any time is 
given by equation (6.2.2). The choice of either <j> or is 
arbitrary as long as this equation is satisfied. The con- 
trol system will be designed to maintain a zero yaw angle. 

Making the assumption of one way coupling for choosing 
control gains, let 

C 4 

N = -k^-kgiJj-kgC eV+Wy)+ | (6.3.1) 

then 

/U ) 2 

i|>+ ( C 3 + b 2 k g ) 4^ + k ^ b 2 ip = - ^gb2 ] ( eV+w^ ) (6.3.2) 


If we require a natural frequency U ' and damping ratio 
5 ^ , a good choice for the gains is, 



O) 



k 5 = (2 5 *“*- c 3 ) ‘ 2 SVz + lr C n r Sq d b < 6 - 3 - 3 ) 

„ 2 C ne Sq d b 
k 6 ■ TWJ - ~ 


Let u, = 1.5 rad sec -1 and r = 0.707 

* <f» 

k 4 = 2.25 I = 2.4 x 1 0 5 N m rad -1 
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kn = 2.121 1+2530 C V 
o z n r 

= (2.2 x 10 5 + 2.5x1 0 3 C V ) N m sec rad" 1 

r 

k c = 337 C V N sec rad' 1 
6 "B 

The second of equations (6.1.6) is of second order in 
<j> and is excited by a linear combination of i|>, and (w y +eV). 
We can choose 

C C 

L = -k 7 +-kg$-k 9 (eV4W y )- ^ ^ (6.3.5) 


The equation for roll angle becomes 

^(Cs+bakgU+bgk?* = (y^- - b 3 k g ) (eV+w y ) 


(6.3.6) 

As before, if w. and s. are the required natural fre- 

9 9 

quency and damping ratio respectively, the gains become 


k 7 

' bj ' *x“* 



k 8 

= FJ ( 25 *V C S ) ' 

25 *Vx + 

2V c t p St| d b 

k 9 

c 6 C * p Sq d b 



VF 3 - 5 - 



wi th 





= 1.5 rad sec" 1 

and = 

0.707 

the control 

gains become. 
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Figure 6.2 The Complete Lateral Control System 







k 7 = 2.25 I = .7xl0 5 N m rad" 1 
k g = 2.121 I x +2 5 30 C £p V 

= (.67x10^ + 2530 C^pV ) N m sec rad -1 (6.3.8) 
kg = 337 C^N sec rad" 1 

Stability derivatives C , C n _, C# and Co. are functions 

n P P P 

of flight parameters. These data are not available for 
Bell Model 266 tilting proprotor aircraft. Therefore the 
gains have not been computed. However, the above equations 
show the order of magnitude of these gains. Once these 
stability derivatives are available the gains can be computed 
as a function of the independent variable, x. 

A block diagram of the control system for lateral con- 
trol is shown in Figure 6.2. 
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APPENDIX A 

Summary of Basic Data for Bell Model 266 
Ti 1 ti ng-Proprotor VTOL Aircraft 


General 

Overall length 

Overall width (proprotors 
turni ng) 

Distance between pylon 
pivot points 

No. of engines 

Shaft horse power (min) 


17.0 m 

26.0 m 

14.3 m 
2 

3435 each (59° day) 
3370 each (92° day) 


Helicopter Mode 


Hoveri ng 

cei 1 i 

ng, out of 

3870 

1 m 

(standard day) 

ground 

effect (design TOW) 

2325 

m 

(95 e 

' day) 


Maximum 

speed 


316 

Km 

hr" 1 

( Twi n 

engine) 

(sea 1 

eve! , 

design TOW) 

259 

Km 

hr" 1 

(Sing! 

e engine) 

High Speed 

Mode 






Maximum 

speed 


654 

Km 

hr* 1 

(Twi n 

engi ne) 

(sea 1 

eve! , 

design TOW) 

470 

Km 

hr’ 1 

(Sing! 

e engine) 

Maximum 

speed 


657 

Km 

hr" 1 

(Twin 

engine) 

(3000 

m, des 

ign TOW) 

505 

km 

hr" 1 

(Singl 

e engine) 

Average 

crui se 

speed 

408 

Km 

hr" 1 

(Sea 1 

evel ) 

(desi g 

n TOW) 

443 

Km 

hr" 1 

(3,000 

m) 




500 

Km 

hr" 1 

(7,500 

m) 

Range 



820 

Km 




(desi g 

n TOW 

and payload) 
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Ming 

Area 

Span 

Aspect ratio 
Airfoil section root 

ti p 

Lift curve slope 

Angle between wing zero 
lift line and fuselage 

Angle of attack at stall 

Aileron-flaps area/side 

Flaps area/side 

Wing loading 


Fuselage 

Length 

Maximum breadth 

(less wheel fairings) 

Area (Planform) 


35.5 m 2 
15.0 m 
6.43 

NACA 64-223, modified 
NACA 64-219, modified 
4.28 rad" 1 
4.9° 

axis 

17° 

3.36 m 2 
3.58 m 2 
358 Kg m" 2 


2 

3 

5.87 m 
108.14 m 2 
10.275 m 2 
. 585 m 
.095 

58.6 Kg m -2 
4.53 


16.1 m 
2.15 m 

36.7 m 2 


Proprotors 

No. of proprotors 

No of blades per proprotor 

Radi us 

Disc area per proprotor 
Blade area per proprotor 
Blade chord 
Solidity 
Disc loading 
Proprotor lock number 
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Horizontal Tail 
Area 
Span 

Aspect ratio 
Lift curve slope 
Airfoil section 
Angle of incidence 

Vertical Tail 

Area (including carry- 
through) 

Span 

Aspect ratio 
Lift curve slope 
Airfoil section root 

ti p 


Weights 

Design take-off weight 
Payl oad 


Weight of rotors, pylons, 
engines etc. (Total) 


9.4 m 2 
6.57 m 
4.61 

4.28 rad" 1 
NACA 64-012 
-0.5° 


12.9 m 2 

5.5 m 
2.34 

2.75 rad" 1 
NACA-64-01 5 
NACA-64-009 


12,700 Kg 
2,990 Kg 
3,175 Kg 
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APPENDIX B 


Summary of the Longitudinal Equations of 
Motion and Weighting Matrices in the Quadratic 
Synthesis Technique 

B.l Nonlinear Equations 
dV 

m = - D^-Dy-Dp-Dj^Q-mg sin y + T cos <t>- H sin $ 
dv 

mV j - +Ly^+ Ly+L ~mg cos y + T sin <|> +H cos $ 


dx 

dt 


V COS y 


^ = V sin Y 


D W ' 7 pS W^R^ C 0y^ a R* 


+C. (a'-Aa) sin 

L W 


L W * 7 pS W V R^ _C D w ^“R 


+ C. (a ' -Aa) COS 
L W 


«f) COS (Aa) 

( Aot ) ]x . 78 + \ pS w V 2 C d (a ,6f )x . 22 

W 

» 6 f ) sin (Aa) 

(Aa)]x.78 + pS^V 2 C^ (a')x.22 


a^ = a - Aa 


ft 

Aa = arcsin[ —■ sin (<f>)] 
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Vn = V 2 +4v 2 +4Vv cos ( <J) ) 

K 

a = 6-Y+.085 
a 1 = a+ . 409 6 f 

<)> = n + 9 - y 

is a solution of 

v[V 2 +v 2 +2Vv cos 4> ] 1 ^ 2 = ' 3 ' c ~ 

p 5 R 

C L (a* ) = 4.28 a' 

4.28 ct^ + [£ -.256f+.4096f-a'] 
[J -.256f+.4096f-a^ + ] 

4.28 a s '"[J + . 2 5 6 f - . 4096f+a 1 ] 
[£ +.25fif-.4096f+a^‘] 

C D (ct',6f) = . 025+ . 1 6f 2 + l . 035a 1 2 

? + 2 

= . 025+ .1 6 f ^ + 1 . 0350^ +.62[a ' 
? 2 

= ,025+.l6f <: + l .035a^~ + .62(a^ 



.31 7 + . 20 5 6 f radians 



o' s = -.31 7 + . 205fif radians 


D 


T 



L 



c l t ( “t ) 


c d t (o t 


) 


pV S T C D t ^°T^ 

P V 2 SyD L (ay) 

= 4.28 a T 
= 1.08 [| - a T ] 

= -1.08 [■£■ +a j] 

= . 025+1. 035a 2 
= .1285+. 62( | a T | - . 31 7) 


l°T I 1 *317 rad 
ctj > .317 rad 

oj < -.317 rad 

|°tj| <_ . 307 rad 
cty > .317 rad 


ay = e-y-,01 

Dp “ ^ pV 2 Sp[ . 05+. 25 s i n ^ ( | o p | ) ] 

L F = \ p V 2 S F [ . 25 sin 2 (ap)cos (ap) sgn(op)] 
ap = e-y 

Nacelle lift and drag are neglected. 
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Thrust is determined as a function of the component of 
the free stream air speed perpendicular to rotor tip path 
plane. In the computer program linear interpolation is 
carried out between points spaced 15 m sec" 1 apart 


Normal Velocity 

Max . Ava i 1 abl e Thrust 

0 m sec" 1 

211 ,000 N 

15 

1 53 ,500 

30 

118,500 

45 

92,500 

60 

71 ,800 

75 

57,000 

90 

47,500 

105 

40,750 

120 

36,250 

1 35 

33,000 


Table B-l Thrust as a Function of Free Stream 
Air Speed. 

The longitudinal inplane rotor force is small under 
most flying conditions and is neglected. 


m 

mass of aircraft 

= 

12,700 Kg 

mg 

weight of aircraft 

= 

123,300 N 

S W 

wing area 

= 

35.5 m 2 

A 

S T 

tail area 

= 

16.25 nr 

S F 

fuselage planform 


36.7 m 2 

S R 

area 

disc area of rotors 
(total) 

— 

216.3 m 2 

P 

d e n s i tv of air 
(at S L ) 

= 

1.267 kg m 
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B.2 Linearized Equations 

The nonlinear equations are linearized about the 
nominal path by numerical differentiation. They are 


dZ 

"37 


= F(x)Z + G ( x ) U 


where 


Z = 


F = 


r$v 

6y 

6h 

W 


3f l 

af l 

n 1 

( 9f i . f iL. 

sin 

w 

8y 

0 I 

Ia7“ + V") cos 

T 

3 f 2 

3f 2 

n 1 

f af 2 „ 2f 2\ 

sin 

w~ 

~ 

0 ! 

\W * T-jcOSy. 

T 

0 

0 

sec 2 y 


0 

0 

0 

0 


0 


/3f : 


sin Y / 3f 2 f . . f l\ 

" f 2 tan Y T~) 


for take-off 


G(x) = 


r 3 f 

77 

3 f / 

i 

3n 

0 


1 


3f 

7T 

3f / 

i 

36 

0 


1 


, U = 


6n 

66 


L 0 
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and for landing 


G(x) 


9f l 

9f l 

9n 

9T 

9 f 2 

9f 2 

9n 

FT 

0 

0 


0 



U = 


<5n 

6T 


and 


dV _ 
dx 


f 


1 



2 


B.3 Weighting Matrices 

The following are the values of the weighting matrices 
used in the quadratic synthesis technique. 

For take-off 


A 


00021 

0 

0 

0 

0 

.04 

0 

0 

0 

0 

o 

.00044 cos y 

0 



0 0 


0 


.14V 0 

f. 3 V > 24.5 m sec" 1 

0 -1 
0 V < 24.5 m sec 
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s 


f 


For landing: 


A = 


B 


-1 


and 


S 


x 


f 


.069 

1.2 

.02 

-.0026 

.2 

170 

2.9 

1.1 

CM 

O 

• 

2.9 

.063 

.013 

.0026 

1.1 

.01 3 

.048 - 


0011 0 

f 

o 

o 

0 .04 

0 0 

0 0 .00044 cos 0 

0 0 

0 0 _ 

.082V 0 


I 

o 

ro 

X 

o 

< 

l 



.36 

.065 

-.0036 

-.04 

.065 

.032 

.0044 

-.06 

-.0036 

.0044 

.0072 

-.057 

-.04 

-.06 

-.057 

240 


108 


